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ANALYSIS OF GRAVITY ASSISTED TRAJECTORIES FOR 
SOLAR-PROBE AND DEEP-SPACE MISSIONS 

By John W. Young 
INTRODUCTION 

As planetary missions become more advanced the energy requirements 
and trip times become increasingly large for vehicles following direct 
transfer orbits from Earth. For example, a solar impact probe requires 
a velocity at the Earth activity sphere (hyperbolic excess velocity) of 
nearly 50 km/sec, while an escape orbit from the solar system would 
require a hyperbolic excess at Earth of 12.4 km/sec. An alternate 
approach to these direct -type transfers made feasible by recent advances 
in orbit determination and vehicle control, is to use the free energy 
attainable by passing through the gravitational field of an intermediate 
planet. This maneuver is commonly referred to as "Planetary Flyby" or 
"Hyperbolic Encounter." For example, proper utilization of the gravity 
field of Jupiter can reduce the required hyperbolic excess velocity for 
solar impact by a factor of three over that for a direct transfer. An 
analogous maneuver can reduce the required velocity for an escape from 
the solar system by nearly one third. 

The general concept of planetary flyby can be stated briefly as 
follows. Consider a vehicle proceeding along a trajectory in helio- 
centric space. For certain orbits the vehicle may approach a passing 
planet such that the planetary gravity field exerts the primary pertur- 
bative force on the vehicle. While the vehicle orbit about the sun is 
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elliptic, the approach and departure with respect to the planet is along 
a hyperbolic path. This hyperbolic encounter can change the energy and 
angular momentum of the vehicle with respect to the sun. For example, 
consider the situation shown in figure 1. The vehicle enters the planet 
activity sphere with a heliocentric velocity, V^ g , and a relative 
velocity with respect to the planet, V^, such that it passes behind the 
planet (direct encounter). After having been deflected by the planet 
gravity field, the vehicle departs on a new orbit with heliocentric 
velocity, V 2 S . Since the velocities of approach (V lp ) and departure 
(V 2p ) with respect to the planet are equal, the encounter results in an 
increase in heliocentric velocity for the vehicle, as is seen in figure 1. 
A reverse type maneuver (retrograde), for which the vehicle passes in 
front of the planet, can result in a decrease in heliocentric velocity 
for the vehicle. Energy, of course, is conserved during a flyby since 
the orbit of the perturbing planet is altered in proportion to the masses 
involved in the encounter. 

Numerous studies relating to planetary flyby have been made in the 
past, as will be outlined in a later section of this report. The objec- 
tive of the present analysis is to supplement these studies by investi- 
gating various aspects of the flyby concept which were generally not 
emphasized in previous studies. An additional aim is to analyze in 
detail various flyby orbits to obtain an understanding of the basic 
mechanics involved. 

Two types of flyby missions will be investigated. One such mission 
is the solar probe, with the objective of either achieving a solar impact 



or of placing the vehicle in close proximity to the solar surface. The 
other type mission studied is the Deep Space mission. The objective of 
this mission is to either escape from the solar system or to achieve an 
aphelion distance greater than that possible without employing the flyby 
concept. 

Results will be presented to illustrate the relative effects of 
various perturbing bodies on significant orbital variables such as 
transfer time, momentum and energy changes at the perturbing planet, and 
launch velocity requirements at Earth. Included will be an analysis of 
the effect on over-all mission performance of varying such flyby param- 
eters as approach and departure direction and velocity and distance of 
closest approach (perigee) at the perturbing body. Consideration will 
be given to determining the sensitivity of specified orbital parameters 
to small variations in certain launch and approach conditions with the 
objective of establishing those quantities most critical to mission 
success. 

In addition to the parametric type results previously outlined, an 
iterative, numerical procedure will be used to obtain a limited number 
of optimal missions. This procedure enables such trajectory variables 
as perihelion and aphelion distances following encounter or orbital 
transfer time to be either minimized or maximized subject to constraints 
on appropriate orbital parameters. Constraints could include launch 
energy and orientation of the transfer orbit at Earth, transfer angle 
and time, or certain approach and departure conditions at the perturbing 
body. These optimal missions are of practical importance for energy or 
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time limited missions and are of academic interest since they give 
insight into the basic concepts involved in planetary flyby. 

Although a complete analysis of the mechanics of planetary flyby 
will be given in a later section, the technique applied in the present 

study can be briefly stated as follows. Since no analytical solution to 

\ 

the three-body problem exists, and since numerical integration of the 
three-body equations of motion is very time consuming, a two-body 
analysis was used. Thus, during a mission it was assumed that the 
vehicle trajectory was influenced by only one body at any particular 
time. Therefore, a typical mission consisted of three distinct segments, 
an initial heliocentric transfer from Earth to the activity sphere of 
the perturbing planet, a swingby of the planet under the influence of 
only the planetary gravity field, and a final heliocentric transfer to 
the desired destination. 
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NOTATION 


a, b, c, e semi-major and semi-minor axes, linear 

eccentricity and eccentricity of tre. sfer orbits 

Cr perigee corridor at encounter planet (eq. (21)) 

p f 

Cp o approach corridor at encounter planet (eq. (20)) 

E eccentric anomaly 

h orbital energy constant 

H auxiliary angle of hyperbolic orbit 

K l' constants 

M mean anomaly 

r radial distance 





time 

orbital velocity 
velocity 

hyperbolic excess velocity of vehicle at 
insertion 

transfer angle associated with vehicle orbit 
elevation angle of velocity vector with respect 
to local horizontal 
true anomaly 

angle between vehicle -planet and vehicle-sun 
radius vectors (fig. 4) 

angle which specifies orientation of insertion 
orbit at earth (fig. 3) 
angles used in determining flyby parameters 
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gravitational field constant 
radius of activity sphere 
angle used in specifying earth and planet 
positions at various mission times 
angle between vehicle-sun and earth or planet-sun 


radius vectors 


Subscripts. - 


A , P aphelion and perihelion conditions 

d direct transfer orbit 

H refers to parameters of post -encounter orbit 

i insertion conditions 

s heliocentric reference 

p planet reference 

E conditions associated with earth 

f refers to parameters of hyperbolic flyby orbit 

1 approach conditions at planet activity sphere 

2 departure conditions at planet activity sphere 
Multiple subscripts . - When dealing with static quantities such as 

velocity, angular orientation, or radial distance, the first subscript 
generally denotes an initial reference position while the second 
subscript specifies the location to which this position is referred. 

For example, consider: 

V, = heliocentric approach velocity at planet activity sphere 
Is 

(1 denotes approach condition, s denotes heliocentric reference) 
Tjps = true anomaly of heliocentric departure orbit at planet activity 


sphere. 
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For quantities such as transfer angles and times , the first and 
second subscripts denote respectively the initial and final reference 
positions. For example consider: 
a il = ^ rans ^ er angle from insertion to encounter 
tii = transfer time from insertion to encounter 
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SURVEY OF LITERATURE 


A comprehensive review of the mechanics of the flyby maneuver, as 
well as problems associated with the practical application of such a 
mission is given in reference 1. An interesting result of reference 1 
is that, based on orbit determination and trajectory control procedures 
available at that time, it was concluded that certain applications of 
the flyby mission were impractical. 

The advent of improved space navigation and guidance techniques has 
largely overcome the difficulties outlined in reference 1 and has lead 
to a renewed interest in flyby type missions. Since this renewed 
interest has resulted in a large volume of literature relating to the 
dynamics of, and requirements for planetary flyby, the present survey 
will be limited to typical studies covering the broad spectrum of avail- 
able literature. In addition to these references relating to the 
dynamics of planetary flyby, literature will also be cited which details 
the mathematics involved in the optimization procedure employed in the 
present analysis. 

References 2 and 3 present results which illustrate the advantages 
of, and requirements associated with the use of the flyby concept as 
applied to specific space missions. Reference 2 gives a general analysis 
of the use of the flyby concept for achieving a multi -planet flyby 
mission during the 1976-1979 'time period. It is shown in reference 2 that 
the Grand Tour mission is feasible within projected state-of-the-art 
capabilities. Reference 3 discusses the requirements associated with 
round trip, Mars stopover missions which employ a Venus swingby technique. 
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The results are compared with those of a direct Mars mission for launch 
opportunities between 1980 and 1999- The results show that, in nearly 
all cases, the flyby significantly reduces the launch requirements at 
Earth. 

Reference k presents an analysis of flyby missions of the type 
considered in the present study. A Lagrange multiplier procedure was 
also used to obtain various optimal missions. While a more complete 
physical model was used in reference 2 than in the present study (planet 
ephemeris as well as out -of -ecliptic transfers were included), the 
results were of a more specific nature since transfers were limited to 
those associated with particular Earth launch opportunities. It also 
appears that the constraints imposed on some transfer orbits lead to an 
omission of certain significant results. 

An analysis of the use of a flyby of Jupiter for achieving Earth 
to Saturn missions is given in reference Results were given to illus- 
trate the increased payloads available using a flyby at Jupiter. 

Reference 6 compares navigational and system capabilities with 
mission requirements for missions such as those of references 2 - 5 as 
well as those of the present analysis. It concludes that relatively 
simple supplements to present navigational systems can achieve the 
accuracies required for flyby missions. 

The numerical optimization procedure used in the present analysis 
is based generaly on the accelerated gradient method for parameter 
optimization presented in reference 7. This reference outlines the basic 
theory involved in the optimization procedure. Additional studies which 


10 


present detailed anaJyses of various aspects of the optimization proce- 
dure are given in references 8, 9 , 10, and 11. 

Implementation of the actual computer program, as evolved from the 
techniques of references 7 - 11, was due largely to the efforts of 
I. L. Johnson (ref. 7) at the Manned Spacecraft Center, NASA. However, 
as no references are currently available in the literature, details 
concerning the computer program and its current application can be 
obtained on request from the author to the present report. 


THEORY 


Before proceeding with a development of the governing equations of 
the present analysis, a general description of the overall planetary 
flyby concept will be given to illustrate the various phases involved 
and to introduce the terminology to be used. A typical flyby mission 
is shown schematically in figure 2. While figure 2 depicts a flyby of a 
planet inside the Earth orbit for which the objective was to achieve a 
reduced perihelion distance, the various phases are similar for other 
missions considered in the study. 

As shown in figure 2, the initial phase involves a heliocentric 
transfer from Earth to the perturbing planet. This transfer, which will 
be referred to as the pre-encounter orbit, is established by specifying 
the energy and orientation of the departure trajectory at the Earth 
activity sphere. This departure from Earth will be referred to as 
insertion. The second phase begins with the entry of the vehicle into 
the activity sphere of the planet. The planetary approach conditions 
are determined from the pre-encounter transfer orbit and from the assumed 
orientation of the vehicle with respect to the planet upon entry into 
its activity sphere. Following a hyperbolic swingby of the planet, the 
final phase begins as the vehicle exits the planet activity sphere. 

This phase, which will be referred to as the post-encounter orbit, 
consists of a heliocentric transfer to the final destination. Orbital 
parameters for the post -encounter phase are determined from exit condi- 
tions at the perturbing planet. 
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The assumptions and limitations of the present analysis as well as 
the physical and mathetical relationships used in the various phases 
shown on figure 2 will now be given. 

Assumptions 

Since the present study is of a general nature as opposed to 
specific studies previously cited in the SURVEY OF LITERATURE, certain 
simplifying assumptions were made. Planet ephemeris was not included 
and planets were assumed to move about the sun in circular orbits. In 
addition the orbit of planets, as well as all transfer maneuvers, were 
assumed to be coplanar. To reduce computational time, a two body analy- 
sis was assumed for all orbital transfers. Thus, a typical flyby mission 
consisted of a series of "patched conics", such that the vehicle was 
influenced by a single gravity field during any particular segment of 
a mission. The assumed physical constants for the planets are given in 
table I. 

Due to the large number of missions available for consideration, 
certain limiting assumptions were also necessary. Only flybys of the 
planets Venus and Jupiter were included. Insertion at Earth was limited 
to directions either along (© 1 = 180°) or opposed (0j_ = 0°) to the Earth 
orbital direction. Thus, at insertion, the vehicle was near the aphelion 
position of the pre-encounter orbit during transfers to Venus. Likewise, 
transfers to Jupiter were initiated near perihelion on the pre -encounter 
orbit. These conditions are near optimal from an energy exchange stand- 
point. One additional limitation on transfer orbits was made to reduce 
the spectrum of possible missions. Pre-encounter orbits were limited 
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to transfer angles of less than 180°. Thus, Venus encounters occurred 
before the perihelion position of the pre-encounter orbit was reached, 
while Jupiter encounters were pre -aphelion. 


Pre -Encounter Transfer Orbit from Earth 
The geometry used to establish the transfer orbit from Earth to the 
perturbing planet is shown on figure 3 . The orbit is determined by 
specifying the hyperbolic excess velocity vector at Earth. Thus, for 
assumed values for velocity magnitude and direction (V he and 0^) at the 
Earth activity sphere, the required conditions for the heliocentric 
orbit can be determined as follows. 

Initial conditions at Earth .- Vehicle distance from the sun at 
Earth is determined by applying the cosine law 


Thus 


r i = 4> + P E ‘ 2r Es P C0S(9 ° ‘ ®i ) 


r i ' < r L + 4 - 2r Es p E Sin * 1 > 


( 1 ) 


Using the law of sines, the angular displacement of the vehicle from the 
Earth-sun line is 

sin = sin ( 90 - 0 ^ 

PZ ri 


or 


♦. = sin -1 (V° S 4 

\ ^ ) 


( 2 ) 


Ik 


The heliocentric 
shown in figure 3> is 


velocity and flight path angle of the vehicle, as 
given by applying the cosine law 



(3) 


and 



00 


Orbital parameters for pre -encounter transfer . - Using the previously 
established initial conditions for the transfer ellipse, the character- 
istics of the orbit can be established. Perihelion and aphelion condi- 
tions can be determined by applying conservation of energy and angular 
momentum. Conservation of energy yields 


V i V P" * 


2 

d 


*i 2 r P d 


conservation of angular momentum gives 


or 


V r cos7 = V p r p 
11 1 r d r d 


cos 


Equations (6) and (5) yield 


( 5 ) 


( 6 ) 


V i- 2 "s V p!’ 2 "s V Pc 
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or 



2u 


V r cosy 




= 0 


Thus, vehicle velocity at perihelion for a direct transfer from Earth is 
given by 




Y ± t ± cob7 ± 


Vicos?, 



(7) 


using this value for V-, , the perihelion distance is found from equa- 

*d 

tion (6). 

In a similar manner aphelion velocity and distance for a direct 
transfer from Earth, can he found to he 



's 


V P r P 
r d 



( 8 ) 



r V 
P d P d 


( 9 ) 


Standard conic relationship (reference 1) for the semi -major and semi 
semi -minor axes and the numerical and linear eccentricities of the 
transfer ellipse are given in equations (A-3) of Appendix A. 
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Hyperbolic Encounter at Planet 

The planetary flyby phase of a mission included the approach to 
encounter at the activity sphere of the perturbing body, the hyperbolic 
orbit about the planet, and the exit from the planet following encounter. 
The equations defining vehicle motion during these portions of the flyby 
will now be developed using the geometry of figure 4. 

Approach conditions at planet activity sphere .- The technique used 
to establish approach conditions at the planetary activity sphere was to 
first establish heliocentric conditions at this point and then to trans- 
fer from a solar to a planetary reference system. The orientation of the 
approach orbit on the activity sphere of the planet, which is needed to 
calculate vehicle radial distance from the sun, was specified by assuming 
a value for the angular displacement (0 ) between the vehicle-sun and 

Cl 

vehicle -planet radius vectors. As will be shown in a later section of 
the report, an assumed value for © a in effect determines the transfer 
time from insertion to encounter and thus establishes the launch date or 
angular separation of the planet and Earth at insertion. 

As shown in figure 4, the angle between the vehicle -sun (r^ g ) and 
planet-sun(r ) radii is given by 


as 


sin 


-1 1 P P Sln 9 a 
r ps 


( 10 ) 


The radial distance of the vehicle from the sun is now given by 


r. 2 = r 2 + p 2 - 2r p cos 180° - + 9 fi ) 

Is ps p ps p Is a 
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or 


Is 




2 2 

r — + P p + 2 r ps p p cos (t 


ps 


Is 


0 ) 

a 


1 1/2 


(ID 


Applying conservation of energy and angular momentum for the pre- 
encounter orbit, the heliocentric velocity and flight path angle at the 
activity sphere are given by 



( 12 ) 


and 




cos 


-1 



cos y ^ 
ls r ls 


( 13 ) 


where is determined by the pre-encounter orbit, as will be explained 

in a later section of the report. 

A coordinate change from a heliocentric to a planetary reference 
can now be made for the vehicle orbit. From figure b it is seen that the 
relative velocity of the vehicle with respect to the planet is 


= t 


IP P 


’2 2 11/2 
U + V, - 2U V. cos(V lB - y. ) 
p Is p Is Is Is' 


W 


also 


0 2 - r lp + 9 a + *ls - 180 ° 


where 


© 2 = sin 


-1 Pis Sln( *ls - 7ls> 
v lp 
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thus 



e + v 


Is 


e 2 - i8o° 


(15) 


Hyperbolic orbit at perturbing planet . - Having determined the 
velocity and flight path angle of the vehicle with respect to the planet, 
the procedures of equations (5) through (9) yield the perigee condition 
at the planet. These are given in equations (A-6) of appendix A. 

Standard conic relationships (ref. l) for the semi-major and semi-minor 
axes, numerical and linear eccentricities, and the true anomaly of the 
vehicle are given in equations (A-7) and A-8) of appendix A. 

Departure conditions at planet activity sphere . - Since the hyper- 
bolic approach and departure conditions with respect to the planet must 
be the same we have 

V 2p = V lp' 7 2p = 7 lp 


In order to determine the orbital parameters for the post -encounter 
transfer orbit, a change must be made from a planet to a solar reference 
system. Referring to figure 4 we have for the vehicle distance to the 
sun 


2 

2s 




- 2p r cos(0, 
p ps '3 



where 

0 3 = 2T1 lf ” 0 1 



and 
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thus 


’2s 


2 2 1 1/2 

= I + r „ + 2p r cos 0 


= P 


ps 


P ps 


( 16 ) 


where 

0 = 2ti + 0 + i|l 

H 'if a Is 

From figure 4 it is seen that the angle between the planet -sun and 
vehicle -sun lines is given by- 


sin sin(9^ - 90°) cos 0^ 


’2s 


2s 


or 


*2s = " sin 


-1 


-JL_ sin 0 
r 2s 


(17) 


Figure 1+ shows that the heliocentric velocity of the vehicle following 
encounter is given by 

V 2s - + 4 ' ^Vp C ° S [ r 2P ' ( ®5 - 90 ° ) J 

or 

V 2S = [ V 4 + ^ + 2V 2p U p COS ( P - (18) 

The heliocentric flight path angle, y , can be found from figure 4 as 


Sin [ 7 2p " ( 0 5 ‘ 9°°J sin ( _ 7 2b - *2s^ cos ( e 3 ‘ >2p) 


2s 


2p 


2s 


follows 
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or 


- sin(7 0 „ + 


2s 


cos(2t) 

2s V If 

2s 


90 + + 

Is 


e - 7 ) 

a 2p 


thus 


= " Sln ’ 1 fe Sln(3 ' W] ' *2s (19) 

Post -Encounter Transfer Orbit 

The post-encounter orbit is, respectively, elliptical, parabolic, 
or hyperbolic as the orbital energy constant, hgg* is negative, zero, cr 
positive. For elliptic orbits the perihelion and aphelion conditions 
can be determined as was previously done in equations (5) through (9) • 

For parabolic or hyperbolic orbits, an escape from the solar system is 
achieved. The conditions which determine the nature of the post- 
encounter orbit are listed in appendix A and the governing equations 
and conic relationships for the orbit are given in equations (A-ll) and 
(A-12). 


Transfer Times and Angles 

Various transfer times and angles of interest for flyby missions 
under consideration will now be established with the aid of figure 5* 
Heliocentric transfer from Earth to planet inside Earth orbit , - 
Transfers to planets inside the Earth orbit were limited to cases for 
which insertion at Earth occurred at or past aphelion (r ) on the 
transfer orbit. Thus, the true anomaly (t^) of the vehicle at insertion 
is less than 180° and the vehicle is initially on the descending leg of 
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the transfer ellipse (negative 7 ). The various transfer angles, times, 

and initial positions of the Earth and planet at insertion and encounter 

are easily obtained from figure 5, and are given in equations (A-13) - 

(A-l 8 ) of appendix A. Note from equations A-k, A-1J, A-16, and A-18 

that by fixing © , the launch date or angular separation of the planet 

and Earth at insertion is established. 

Heliocentric transfer from Earth to planet outside Earth orbit . - 

Transfers to planets outside the Earth orbit were limited to cases for 

which insertion at Earth occurred at or past perihelion (r_ ) on the 

a 

transfer orbit, as is shown in figure 5. Thus, the vehicle at insertion 
was on the ascending portion of the transfer ellipse (+7 ). The 

geometry involved for these types of transfers is given in figure 5 
and the various relationships of interest are listed in equations (A-13) - 
(A -18) of appendix A. 

Post -encounter transfer orbit .- As was previously stated, the post- 
encounter orbit is dependent on the orbital energy constant h 2g , being 
elliptic, parabolic, or hyperbolic as h 2 s is, respectively, negative, 
zero, or positive. Since an escape from the solar system is achieved 
if the post -encounter transfer is hyperbolic or parabolic, no transfer 
angles or times were computed for these orbits. 

Two possibilities exist for elliptic post -encounter orbits since the 

heliocentric flight path angle ( 7 . ) can be either positive or negative. 

2 s 

The geometry used to establish the transfer angles and times for this 
phase of a mission is given in figure 6 . Transfer angles and times for 
the post -encounter orbit are given in equations (A-22) and (A-23) of 


appendix A. 
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Transfer time from insertion at Earth to the perihelion or aphelion 
position on the post -encounter orbit is obtained by summing the individual 
times for each phase of the mission and is given in equations (A-21+) of 
appendix A. 

Flyby maneuvers result in a displacement of the major axis of the 
post -encounter transfer orbit from that of the direct transfer orbit. 

It is of interest to establish this displacement in order to determine 
the positions of the Earth and planet with respect to post -encounter 
perihelion and aphelion. The geometry used in defining this is given on 
figure 7 and the desired equations are given in (A-25) of appendix A. 

The equations developed in preceding sections of this report and 
the required standard conic relationships are summarized in appendix A. 
While these equations define vehicle motion from insertion at Earth to 
aphelion or perihelion on the post -encounter orbit, certain constraints 
must be applied to the pre -encounter orbit in order to achieve a desired 
encounter. These conditions will now be outlined. 


Necessary Conditions at Insertion and Encounter 

Insertion conditions . - An obvious constraint on insertion conditions 

is that the pre -encounter orbit must intersect the orbit of the planet 

with which an encounter is desired. Thus, for encounters with a planet 

inside the Earth orbit, the perihelion (r p ) of the transfer orbit must 

*d 

be less than the planet-sun radial distance (r ). Likewise, for 

ps 

encounters with planets outside the Earth orbit, the aphelion (r ) of 

*A 

the transfer orbit must be greater than r^ g . Hence, the hyperbolic 
excess velocity vector (V he ) at Earth must be such that the preceding 
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conditions are fulfilled. These conditions are given in (A-26) of 
appendix A. 

Encounter conditions . - Constraints on encounter conditions at the 

activity sphere of the perturbing planet will be defined with the aid of 

figure 8 . This figure shows how the vehicle is postioned with respect 

to the perturbing planet by specifying the angle (9 ) between the vehicle- 

& 

planet and vehicle-sun lines. Thus, for values of 9 between 0 and 

a 

and 180°, the approach is in quadrant I (0 < 9_ < 90°) or quadrant 
II (90° <9 < 180°) of figure 8 ; while for larger values of 9 the 
approach is in quadrant IH (l80° < Q & < 270°) or quadrant 
IV (270° < 9 a < 360°). 

From figure 8 it is seen that for an encounter to occur, the 

heliocentric conditions for the vehicle must be such that an entry is 

made into the activity sphere of the planet: Thus, the relative velocity 

of the vehicle with respect to the planet (V ) must be directed into 

IP 

the activity sphere (- 180 <7 < 0). These encounter conditions are 

lp 

listed in (A-27) of appendix A. 

It should be noted that, if the flight path angle of the flyby orbit 

( 7 ^ ) is equal to -90°, an impact with the planet will result. Thus, 

the magnitude of 7 , determines whether the encounter is retrograde or 

IP 

direct and the limiting values for each quadrant is apparent from 

figure 8 . For example, a retrograde maneuver occurs for quadrants I 

and IV if, 7 is between 0 and -90° while a direct encounter would 
IP 

result in quadrants II and III for the same approach conditions. 
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RESULTS AMD DISCUSSION 

The results to be presented were, as previously stated, obtained by 
two methods. One method consisted of parametrically varying launch 
energy and orientation (©^, V^ e ) and approach angle (0 & ) at the activity 
sph?re of the encounter planet. The second method of analysis made use 
of a numerical optimization procedure to obtain values for the launch 
and approach variables such that some orbital function was optimized 
subject to constraints on other orbital parameters. 

It should be stated that, while the results could have been obtained 
without reliance on the optimization computer program, its' use greatly 
reduced the data collection task. For example, a specific flyby mission 
satisfying certain criterion can be obtained with a single computer run 
using the optimization program. This same mission might otherwise 
require a large number of runs with a resultant lengthy data reduction 
effort needed to obtain the desired flyby result. Thus, the optimization 
program was, in general, used to reduce computer and data reduction time 
by establishing nominal values for mission variables. The variables 
were then varied parametrically about those nominal values. 

In the following presentation of results no distinction will be 
made regarding the method of data collection. 

Direct Transfer Orbits from Earth 

Direct transfer orbits from Earth were determined to provide an 
index for comparison of the performance of planetary flyby missions as 
well as to establish possible pre-encounter orbits. Parameters of 
interest for these direct transfers are given in figure 9* 
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Shown on figure 9(a) is the perihelion distance and time from 
insertion to perihelion for various values of hyperbolic excess velocity 
at Earth. As previously stated, these transfers inside the Earth orbit 
were established by inserting the vehicle in a direction opposite to the 
Earth's orbit (0^ = 0). Note that insertion velocities of at least 
2.5 km/sec are required for encounters at Venus. 

Figure 9(b) gives the aphelion distance and transfer time for 
launches outside the orbit of the Earth. For these transfers, the 
insertion angle (0^ was 180°. Note that launch velocities of nearly 
9 km/sec are required to effect an encounter at Jupiter. Due to the 
strong influence of launch energy for these transfers, a wide variation 
in transfer time to aphelion occurs, as is seen on figure 9(b). 

Solar Probe Missions 

Results were obtained for solar probe type missions which involved 
intermediate flybys of Venus and Jupiter. The objectives were to 
determine the influence on overall mission performance of various 
launch (V^, 0^) and planetary approach (0 & ) conditions. An additional 
aim was to establish areas of interest with respect to mission optimi- 
zation and to define criterion for evaluation of these optimal type 
missions. 

Venus flyby . - The basic nature of a flyby orbit at Venus will be 
illustrated with the use of figures 10 and 11. Shown on figure 10 is the 
post -encounter, perihelion distance, and the transfer time from insertion 
to perihelion for various launch and Venus approach conditions. The 
vehicle was inserted at Earth in a direction opposite the Earth's orbital 
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direction (0.^ = 0), and the encounter condition was pre -per ihelion 
(see fig. 5). As is seen on figure 10.; perihelion and transfer time 
are strongly influenced by the approach condition. The discontinuity 
in the curves represents conditions for which an impact at Venus occurs. 

Consider the variation of perihelion distance with 0 given on 
figure 10 for V he = 6 km/sec. As 0 & is increased from zero, peri- 
helion increases from the direct transfer value of about 0.468 AU to 
about 0.6 AU for 0 D = 9*5°. For this range of approach conditions, the 
encounter is direct (see fig. l). The increase in r^ with 0 & is a 
result of the vehicle approaching closer and closer to the planet surface 
as 0 approaches its limiting value of 9*5°. This limiting value 
represents the condition for which the vehicle passes the planet at a 
distance of one planetary radii (6200 km) . While the one planetary 


radii limit was chosen for computational convenience, in an actual 

mission the approa*: would be limited to about 1.1 planetary radii so 

as to remain above the planet atmosphere and thus prevent possible 

contamination or impact with the planet. 

As 0 fi is further increased, the encounter orbit passes from one 

side of the planet to the other, and at 0^ = 11° the distance of 

closest approach is again one planet radii. For this condition, the 

encounter is retrograde and a minimum value for perihelion is achieved. 

For large values of 0 , flyby distance increases and r p approaches 

a r H 

the direct encounter condition. The upper limit for 0 of about 85° 

a 

shown on figure 10 corresponds to the condition for which the vehicle 
just intersects the activity sphere of the planet. For larger values 
of 0 a , encounters are not possible. 
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As is shown on figure 10, an increased launch velocity yields a 
reduced perihelion. The transfer time curves shown on figure 10 illus- 
trate a basic fact associated with most physical problems. That is, 
increased performance with respect to one variable must generally come 
at the expense of some other variable. For example, as seen on fig- 
ure 10, a reduced perihelion yields an increased transfer time and 
vice versa. 

Consider again the two limiting cases for direct (0_ = 9-5°) and 

retrograde (0 a = 11°) encounters shown on figure 10. To obtain a 

physical picture of the mechanics involved, it is instructive to examine 

the actual flyby orbits for the two cases. This is illustrated 

schematically on figure 12 for the two limiting encounters. Given are 

the approach and departure conditions at the planet activity sphere, as 

well as additional flyby parameters. While radius vectors given on 

figure 11 cannot be shown to scale, velocities and angles are properly 

scaled. Note that the encounter conditions for both cases are almost 

the same, the only significant difference being that elevation angle 

(y ) for the direct or retrograde cases is respectively slightly 
IP 

greater than or less than 90°. Thus, in the direct encounter the vehicle 
passes behind the planet while for the retrograde case the vehicle 
passes in front of Venus. 

While approach angle to the planet (0 & ) is a convenient variable 
mathematically for establishing the encounter orbit, a more practical 
parameter for evaluating the orbit is the distance of closest approach 
(r ) to the planet. Shown on figure 12 is the solar perigee attainable 
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as a function of distance of closest approach at Venus for various launch 
velocities at Earth. For comparison purposes, the perihelion distance 
for direct transfers from Earth is also given. From this figure it is 
seen that as the distance of closest approach at Venus is reduced, the 
post -encounter perihelion is also reduced, with a maximum decrease of 
about 65 percent (over that for the direct launch) occurring at a flyby- 
perigee of one planetary radius. 

Total mission time to perihelion and Earth position with respect 
to perihelion are given on figure 13 for the range of conditions included 
on figure 12. For comparison, direct transfer conditions are also 
included. As shown, the flyby maneuver only slightly effects total 
mission time. From a pratical standpoint, the increased times required 
for the flyby orbits would have little effect on overall mission 
performance. 

The variation in Earth position with respect to post -encounter 

perihelion given on figure 13, is of practical importance since the Earth 

should be in a favorable position to receive telemetry transmission from 

the vehicle during perihelion passage. While the ideal condition might 

be for the Earth and vehicle to be aligned (0-gpjj = 0), it is seen on 

figure 13 that line of sight contact is achieved for the conditions 

considered and that the increase in transmission distance is not 

excessive. For example, consider the values given for =6 km/sec 

he 

and r_ = 1 planet radii. The increased transmission distance for the 

p f 

flyby mission is about 50 percent greater than that for the ideal case 
(^EFH = 8X1(1 only at) ° ut ^ P^cent greater than for the direct transfer 


orbit. 
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Jupiter flyby .- In the preceding analysis of Venus flyby, solar- 
probe missions, the results presented followed a somewhat predictable 
pattern. However, since the gravitational field at Jupiter is much 
larger than that at Venus, solar-probe missions involving flybys at 
Jupiter are more interesting and unpredictable in nature, as will be 
shown in the following discussion. 

Figure 14 shows the variation in solar perihelion with distance of 
closest approach at Jupiter for several values of launch velocity. The 
results of figure 14 were for a launch along the Earth orbital velocity 
vector (0-j_ = 180°). In contrast to the results for similar flybys at 
Venus (fig. 12), minimum perihelion is not attained at the limiting 
flyby distance of one planetary radius and, for some launch velocities, 
multiple minima occur. Also, while the minimum perihelion for a Venus 
flyby was about 0.146 AU (V^e = 12 km/sec), solar impacts are actually 
attained for the minima shown on figure 14 associated with launch 
velocities of 11 and 12 km/sec. In order to explain the existence of 
multiple minima shown on figure 14, it is of interest to examine the 
basic mechanics of the flyby orbits involved. Consider the minima 
attained for V he = 11 km/sec at distances of closest approach of J.9 
and 12.9 planetary radii. The encounter orbits are depicted on fig- 
ure 15. The velocities and angles shown on figure 15 are to scale while 
the radial distances are not. 

Note that encounter conditions are nearly the same for each orbits 
shown on figure 15 , the only significant difference being that the 
initial flight path angle is about one degree nearer the vertical for 
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the closer flyby orbit. This closer approach gives a greater deflection 

angle at the planet with the result that, following the encounter, the 

vehicle's heliocentric velocity vector is practically aligned with the 

radius vector to the sun (7 2s 63 -90°). For the r p ^ = 12.9 case, the 

smaller deflection angle results in a positive heliocentric flight path 

angle of about 90° as shown on figure 15. Thus, the vehicle proceeds 

outward along the vehicle-sun line and then falls back into the sun. 

A similar effect to that shown on figure 15 accounts for the minima 

shown on figure 14 for V he = 12 km/sec. Note that for this case there 

are actually three values of r^ for which minima exist, with two values 

P f 

occurring for r <2. These lower minima result as the post-encounter 

P f 

flight path angle (y ) passes successively through -90 u in an increasing 

and decreasing fashion. While the scale on figure 14 does not show it, 

a similar effect occurs on the V he = 11 km/sec curve between values 

for r of 5 and 4 planet radii. 
p f 

The single minima shown on figure 14 for launch velocities of 9 
and 10 km/sec is due to a different effect from those previously 
described. These result not from the vehicle being directed either 
toward or away from the sun, but from a reduction in heliocentric 
velocity such as occurred for retrograde encounters at Venus (fig. 11). 

A further examination of figure 14 shows that, for = 9 km/sec, 

the direct transfer condition yields a lower perihelion than on the flyby 
orbit. Thus, only velocities of at least 10 km/sec will be included 
in further analysis of Jupiter, solar-probe missions. 

Based on perihelion alone, the Jupiter flyby is far superior to the 
Venus flyby, as a comparison of figures 14 and 12 show. Also, for the 
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Jupiter case, solar impacts can be achieved at greater flyby distances, 

thus reducing the risk of planetary impact. Again, considering only 

perihelion, the minima corresponding to flybys at about 11 planetary 

radii are preferable to those at closer flyby distances. However, a 

complete evaluation of the Jupiter, solar-probe mission must include an 

analysis of other significant trajectory variables. 

An important factor for flyby missions is the total trip time. 

Mission times for the encounter conditions discussed in figure 14 are 

given on figure 1 6. This figure shows that, for the launch velocities 

included, time from insertion to perihelion increases almost linearly 

with distance of closest approach at Jupiter. Thus, for V =12 km/sec, 

he 

mission time increases from about 950 to 2900 days respectively for the 

flyby distances of 1.5 and 14. 5 planetary radii corresponding to minimum 

perihelion conditions. Mission times of about 1050 and 1700 days occur 

for the lower and upper minima in the V, =11 km/sec case. This 

he 

increase in mission time with flyby distance is in general, caused by 

the effect described on figure 15. That is, as r becomes greater, 

P f 

the post -encounter flight path angle (? 2s ) increases from negative to 
positive values such that the post -encounter orbit progresses from a 
pre -perihelion to a pre-aphelion condition. 

Consider now the position of the Earth with respect to perihelion 
passage by the vehicle. This quantity, which is of practical signifi- 
cance from a telemetry standpoint, is shown on figure 17 for the Jupiter 
flyby missions in question. Since this displacement angle increases 
nearly linearly with time, for simplicity only the magnitude of the 
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angle and those values near minimum perihelion are included on figure 17. 
Positive angle displacements from perhelion are found on the negative 
slope portion of the curves on figure 17. 

Consider the = 12 km/sec data of figure 17. Note that for 

the close planet flyby situation, the Earth and vehicle are at near 
opposition during minimum perihelion conditions. Thus, telemetry trans- 
mission would be impractical in this region. The close flyby, minimum 
perihelion conditions for V he = 11 km/sec are seen to yield very 
vavorable telemetry conditions. To a somewhat lesser degree, this is 
also true for all minimum perihelion orbits at the larger flyby 

distances (r > 11 planet radii). 

P f 

An analysis of the Jupiter flyby, solar-probe mission based on the 
^ata of figures 14, 1 6 , and 17; that is, distance of closest approach to 
the planet, mission time, and Earth position at perihelion passage, leads 
to the following conclusions. As previously stated, from a consideration 
of reduced risk of impacting the planet, flyby distances greater than 
about 11 planetary radii are desirable (fig. 14). Telemetry transmission 
conditions are also favorable for these flyby distances (fig. 17). 
However, the increased times required for these missions (fig. 1 6), 
particularly for the highest energy launch, may necessitate flybys at 

the minimum perihelion conditions achieved for r <5 planet radii. 

P f 

Under these circumstances, the lower energy launch (V hg = 11 km/sec) is 
preferable since a solar impact is achieved at a greater flyby distance 
in addition to having the Earth in a near ideal position for telemetry 
transmission. It is also of interest to note that, for the 
V he = 11 km/sec launch, the perihelion curve is nearly flat for planet 


33 


approach distances between about 3 and 15 planet radii. Therefore, this 
transfer orbit would appear to be less sensitive than the others to 
errors in approach conditions. However, a final evaluation must include 
an error sensitivity analysis for the various flyby orbits. This analysis 
will be given in a later section of the report. 

Deep Space Missions 

In the present study, deep space type missions of a specific nature 
were not considered. The analysis was limited to establishing the maxi- 
mum performance attainable for such missions through use of an inter- 
mediate flyby of Venus and Jupiter. 

Venus flyby . - Various flyby orbits at Venus were computed to 
determine the max imum post -encounter aphelion attainable. The results, 
which were in agreement with reference 4, showed only marginal improve- 
ments over those for direct transfers from Earth. Thus, no data relating 
to Venus flyby - deep space missions is included in the present report. 

Jupiter flyby .- In contrast to Venus flybys, encounters with Jupiter 

provide a definite advantage over a direct Earth launch for deep-space 

missions. This capability is shown in figure 18. This figure gives the 

* 

maximum heliocentric velocity attainable following a direct encounter 
as a function of distance of closest approach at Jupiter and excess 
velocity at Earth. The heliocentric velocity which would result at the 
orbit of Jupiter on a direct transfer from Earth is also given for 
purposes of comparison. 

Figure 18 shows that, for a wide range of flyby orbits, significant 
increases in heliocentric velocity are attainable. An escape from the 


solar system is, in facx, achieved for the results shown on figure 18, 
since the solar system escape velocity at Jupiter is about 18 km/sec. 

The required direct launch velocity (v^ e ) at Earth to effect a solar 
system escape would be about 12. k km/sec. 

While detailed results relating to flybys at Jupiter for use in 
other outer solar system missions will not be given, one such comparison 
is of interest. Consider the case of a Saturn mission. A direct launch 
to Saturn would require a hyperbolic excess velocity at Earth of at 
least 10.1+ km/sec. For this launch, the heliocentric velocity at 
Jupiters orbital distance would be about 13.5 km/sec. Thus, it is seen 
from figure 18 that a Saturn mission could be achieved with a launch 
velocity as low as 9 km/sec by using a Jupiter flyby. 

A feature of interest on figure 18 is that the maximum post- 
encounter velocity occurs at a flyby distance of about 2 planet radii. 
This results since the deflection angle at Jupiter is such that the 
vehicle exits in a dfrection nearly aligned with the orbital velocity 
of Jupiter (7 2g ~ 0). Thus, for this flyby distance, a maximum increase 
in heliocentric velocity is realized (see fig. 11). 

Guidance Analysis for Flyby Missions 

Various questions of practical significance arise in the analysis 
of results presented in previous sections of this report. An evaluation 
of the idealized data given may lead to premature conclusions regarding 
the relative merits of a particular mission. Thus, any final analysis 
must consider the guidance accuracy requirements as compared to that 
attainable with current state of the art techniques. In addition, the 






sensitivity of critical orbital quantities to variations in launch and 
approach parameters must be included. A complete sensitivity study of 
the flyby concept as applied to particular mission applications would 
require a thorough error analysis and is beyond the scope of the present 
study. However, certain conclusions of a general nature regarding the 
guidance requirements for flyby missions can be made. 

Approach corridor analysis . - A standard procedure for establishing 
the guidance requirements of a particular space mission is to specify a 
corridor width at some position on the transfer orbit. This corridor 
represents a "keyhole" through which the vehicle must pass to accomplish 
the desired mission. In any detailed analysis, the corridor at some 
point in space can be pictured as an "error ellipsoid" . The axes of the 
ellipsoid would specify allowable deviations in position and velocity 
components. Since the present study is limited in scope, a simplified 
corridor analysis will be used. 

Two accepted procedures for specifying planetary approach corridors 
are to give allowable deviations in either the vehicle position at some 
distance from the planet, or in the perigee distance achieved during 
passage at the planet. Both methods will be used in the present study. 

Allowable deviations in position will be determined at the activity 
sphere of the planet. This corridor, which represents a segment of arc 
which must be intersected at the activity sphere by the vehicle orbit, 
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where 50 a is a small increment in the approach position on the activity- 
sphere, and is the resulting change in vehicle -sun-planet angular 

displacement. The perigee corridor is merely the difference in perigee 
attained for the two approach conditions shown on figure 19 and is given 


by: 



( 21 ) 


where the unprimed and primed quantities refer respectively to the 0 

& 

and 0 + S 0 D approach conditions. 

cl a 

This simplified representation neglects certain factors which would 

have to be considered in a more complete analysis. While the small 

increment in 0 yields nearly identical flyby orbits, there are minute 

differences in the approach velocity and elevation angle (V ,7 ) . 

Ip lp 

Although inclusion of these small differences in approach conditions at 

the corridor limits would slightly alter the resulting corridors, the 

changes are insignificant as applied to the present study. 

Venus encounter corridors . - Shown on figure 20 are variations in 

approach and perigee corridors at Venus with percent of minimum solar 

perihelion attained. Zero corridor conditions on this figure represent 

the case for which the vehicle passes Venus at a flyby distance of one 

planet radius (fig. 12). Thus, as the corridor increases, perihelion 

and flyby distance (r ) also increase. As seen on the figure, launch 

pf 

velocity at Earth has little effect on approach and perigee corridors at 
Venus. 

From a practical consideration, figure 20 can be interpreted as 
follows. Assume that, for a particular solar probe mission; launch, 
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orbit determination, and guidance errors are such that the vehicle can 

be positioned on the activity sphere of Venus within a 4000 km distance 

of the aim point. Thus, to achieve minimum perihelion while avoiding an 

impact at Venus, the approach corridor (Cp ) would be 8000 km with the 

a 

aim point being the 4000 km corridor condition. Likewise, the perigee 
corridor would be one planet radius in width with the aim point being at 
a flyby distance of 1.5 planet radii. For this situation, solar 
perihelion would be between 60 and 100 percent of the minimum achievable 
value. 

As was shown in reference 6, expected accuracy in delivering the 
mariner-Venus probe was within a 100 x 250 km ellipse at the target aim 
point. This accuracy was to be achieved with only Earth based tracking 
and did not rely on on-board measurements which could be expected to 
further improve planet approach guidance. Thus, assuming similar 
approach guidance accuracy, figure 20 shows that a successful Venus 
flyby-solar probe mission could be achieved with at least 95 percent of 
minimum perihelion being attained. 

Jupiter encounter corridors .- Figure 21 gives the variation in 
approach and perigee corridors with achievable solar perihelion for 
various launch velocities at Earth. The results shown are for conditions 
near the multiple minimum perihelion values occurring at flyby distances 
less than and greater than 6 planet radii (fig. 14). 

Corridors associated with the larger flyby distances are given on 
figure 21(a). As before, zero corridors occur at the minimum perihelion 
condition shown on figure lb. Positively and negatively Increasing 
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corridors refer respectively to flyby distances greater and less than 
the minimum perihelion flyby distance. 

An analysis of figure 21(a) shows that the corridors decrease as 
launch velocity is increased. Positive corridors are seen to be some- 
what larger than corresponding negative corridors. In contrast to Venus 
flybys (fig. 20) for which the aim point cannot be the zero corridor 
case, the optimal aim point for a Jupiter flyby would be at the zero 
corridor condition. 

Shown on figure 21(b) are corridors for the lower flyby distance - 
minimum perihelion missions shown on figure 14. While previous figures 
have shown zero corridors at the 100 percent condition, figure 21(b) 
reveals a sizeable corridor width at this point. This results since the 
perihelion curves of figure 14 are fairly flat at the lower minimum flyby 
conditions and a solar impact is actually attained over a range or orbits 
about the minimum. 

As seen on figure 21(b), the lower energy launch gives significantly 

larger approach and perigee corridors than for the higher energy launch. 

Also, the chance of planetary impact is greatly reduced for the 

V =11 km/sec transfer since minimum perihelion is achieved at a 
he 

greater flyby distance. Thus, consideration of flyby corridors, transfer 
times (fig. 1 6), and telemetry transmission conditions (fig. 17) reveals 
a distinct advantage for the lower (V he = 11 km/sec) as opposed th the 
higher (V he = 12 km/sec) energy transfer orbit. 

In order to establish approach and perigee corridor conditions for 
which a solar impact is achieved, various flyby orbits were computed for 
launch velocities between 10 and 12 km/sec. The results are shown on 
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figure 22. The corridors given are equivalent to the IOC percent corri- 
dors shown on figure 21(b). Note that maximum corridor width occurs at 
a launch velocity of about 10.8 km/sec and that, for launch velocities 
as low as 10.5 km/sec solar impact is achieved with wider corridors than 
those for launch velocities greater than about 11 km/sec. 

The variation of certain flyby parameters for the maximum corridor 
condition of figure 22 is of interest and is given on figure 23 . Shown 
is solar perihelion distance, transfer time, post -encounter heliocentric 
flight path angle and velocity, and vehicle -Earth angular separation at 
perihelion passage for various flyby distances at Jupiter. An analysis 
of figure 23 shews why a solar impact is attained over such a wide range 

of flyby distances. For r >8.5 and r <6.5 planet radii, post- 

f f _ 

encounter flight path angles are close to + 90 ° and - 90 ° respectively. 

Thus, small perihelion values are attained. For encounter conditions 

such that 8.5 > r-j >6.5 the reduced perihelion results from low values 

for post -encounter, heliocentric velocity. This variation in V 2g and 

7 results since the vehicle velocity vector with respect to Jupiter 

(V ) is nearly opposed to the planet velocity vector following encounter 
2p 

(see fig. 15 ). 

It is apparent from previously given results that, with respet to 
overall mission performance, many tradeoffs exist between such flyby 
parameters as launch velocity, mission duration, encounter corridors, 
and telemetry transmission conditions. It is of interest to note that, 
for Jupiter flyby-solar probe missions, superior overall mission per- 
formance is attained with launch velocities less than the maximum 


assumed values. 
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The encounter corridor analysis has been limited to solar probe 
missions and has not included deep space type missions. No corridor 
analysis for deep space missions is included since the study of these 
transfers was limited to solar system escape oibits. Also, for this 
type mission, post-encounter conditions are fairly insensitive to flyby 
conditions, as was shown on figure 18. Had consideration been given to 
specific deep space missions than a guidance corridor analysis would be 
required. For example, a Jupiter flyby-Saturn mission with an objective 
of placing the vehicle between the rings of Saturn puts rather stringent 
limits on approach guidance. 


SUMMARY AND CONCLUSIONS 


A study has been made of the planetary flyby concept as applied to 
solar probe and deep space missions. The objectives were to supplement 
previous studies and to analyze various aspects of the flyby mission to 
acquire an understanding of the basic mechanics involved. 

A planar two-body analysis was used in the study. Thus a typical 
mission consisted of a pre -encounter, heliocentric transfer from Earth 
to the activity sphere of the perturbing planet, a planet referenced 
flyby orbit, and a heliocentric transfer the final destination. An 
iterative, numerical procedure was used to obtain various optimal 
missions subject to constraints on appropriate orbital parameters. 
Variables included in the study were the energy and orientation of the 
launch trajectory at Earth and the approach direction to the perturbing 
planet. Hyperbolic excess velocities at Earth between 6 and 12 km/sec 
were used in the study. 

Results are presented to illustrate the effect of flybys at Venus 
and Jupiter on such parameters as transfer time, momentum and energy 
changes at the perturbing planet, approach and flyby corridor widths, 
and conditions at perihelion passage with respect to possible telemetry 
transmission from the vehicle to Earth. 

The results of the study can be summarized as follows: 

1. Maximum reductions in perihelion of about 65 percent over that 
for direct transfers are possible for Venus flyby-solar probe missions. 
Only small increases in mission times result. Favorable vehicle to 
Earth telemetry conditions exist at perihelion passage for these missions. 
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Approach and flyby corridors at Venus were found to be such that, based 
on current guidance techniques, at least 95 percent of the minimum 
perihelion achievable under ideal conditions can be realized. 

2. For Juipter flyby-solar probe missions with launch velocities 
greater than about 10.5 km/sec, solar impacts are possible for multiple 
flyby distances between 1 and 15 planet radii. Numerous tradeoffs exist 
at these minima between transfer time, telemetry conditions, launch 
velocity, and approach and planet perigee corridors. Superior overall 
mission performance was attained for a launch velocity of about 
10.8 km/sec for the conditions considered in the analysis. 

5. Only marginal marginal improvements over that for direct trans- 
fer orbits were achieved for Venus flyby-deep space type missions. 

4. Significant increases in performance were achieved for Jupiter 
flyby-deep space missions. Solar system escape orbits are achievable 
over a wide range of approach and flyby distances for hyperbolic excess 
launch velocities as low as 9 km/sec. 
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APPENDIX A 

SUMMARY OF EQUATIONS USED IN ANALYSES OF PLANETARY FLYBY 

Equations used in the analysis are summarized below for each phase 
of the flyby mission. Equations unique to this analysis are derived in 
the Theory section of the report, while standard conic relationships 
used are found in reference 1. 

Initial values are assumed for the variables 0., V, , and 0_. 

1 he a 

Values for various constants appearing in the following equations are 
listed in tables I and II. 


Pre-Encounter Transfer Orbit from Earth 
Heliocentric Departure Conditions at Earth Activity Sphere . - 
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Hyperbolic Encounter at Planet 
Heliocentric approach conditions at planet activity sphere. - 
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Planet referenced approach conditions at activity sphere. - 
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Perigee conditions at planet. 
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Conic relationships for encounter orbit. - 
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True anomaly of encounter orbit at activity sphere . - 

I + 7 ip -iR b f ■ p p a f* 

l§ + 7 ipl [ Vf 

Planet referenced departure conditions at activity sphere. 


V 2p - V 1P 


(A-8) 


(A-9) 


7 2p ' 7 lp 


Post -Encounter Transfer Orbit 

Heliocentric departure conditions at planet activity sphere. - 
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Conditions on post -encounter orbit. - 



h 2 s < Elliptic Orbit 

^2s = ® Parabolic Orbit 

h > Hyperbolic Orbit 







Transfer Times and Angles 


Vehicle true anomaly at insertion and encounter. - 
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Vehicle eccentric anomaly at insertion and encounter. - 
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Vehicle transfer angle from insertion to encounter. - 
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Vehicle transfer time from insertion to encounter . - 
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Earth position with respect to perihelion (r p ) at insertion and 


encounter. - 
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Planet position with inspect to perihelion (r ) at encounter and 

x d 


encounter. - 
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Auxiliary angle and mean anomaly of hyperbolic orbit at planet. - 
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Transfer time in activity sphere of perturbing planet. - 
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Vehicle true and eccentric anomaly following encounter. - 
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Transfer angles from encounter to perihelion and aphelion. - 
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Transfer times from encounter to perihelion and aphelion. - 
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Transfer time from insertion to perihelion and aphelion for direct 


and flyby orbits. - 
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Earth and planet position with respect to post -encounter perihelion 


and aphelion. - 
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? Necessary Conditions at Insertion and Encounter 

Encounter with planet inside or outside earth orbit. - 
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DESCRIPTION OF OPTIMIZATION PROCEDURE 

In the section on SURVEY OF LITERATURE various references are given 
which provide a detailed description of the theory associated with the 
numerical optimization procedure used in the present study. Thus, only 
an outline of the basic procedure and the practical aspects of its 
application will be presented. 

The optimization procedure provides a method for determining the 
values of certain variables x which minimize some function of the 
variables f(x) while satisfying constraints g(x) = 0, where 
x is an n-vector, f is a known scalar function 
of x, and g is an m-vector of known functions 

of x. The general form of the function to be optimized is given by 

f = f + i I G. g. (B-l) 

J = 1 J J 

where the G’ s are weighting factors associated with the constraints and 
are chosen such that a proper balance is maintained between the 
individual contributions to f of the unconstrained function f and the 

constraints g . 

J 

Constraints may be of either an equality or inequality type. In the 
case of an equality constraint, the value of some parameter is required 
to be identically equal to a desired value. For inequality constraints, 
an upper and lower bound is placed on the parameter and it is allowed to 
assume values within these bounds. The notational form of equality and 
inequality constraints are as follows: 
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Equality - g = (gy - g D ) 

Inequality - g = 0 for (lower limit) <_&„.< (upper limit) 

= (g v - lower limit) for g v < lower limit 
= (upper limit -g^.) for g v > upper limit 
where g v is the computed value of the constrained parameter and g^ 
is the desired value. 


It should be noted that, since the optimal value for f is achieved 

in an iteia ve fashion, violations of the constraints may occur during 

certain iteration cycles. However, as the optimum is approached, the 

constraints are forced within the required tolerances. 

In addition to the unconstrained function f and the constraints g, 

the optimization procedure requires certain other quantities in order to 

establish extremum values for f. These include the effects on the 

unconstrained function f and the constraints g. of variations in 

d 

the variables x. The sensitivity of f to changes in x is the 
gradient of f and is given in the form of an array 


/df 


(B-2) 


The effects on the constraints of variations in the variables are 
supplied as an mxn matrix as follows 

Ai \ s g A 

/s*s£-s-\ 


S Sx 2 


(B-3) 
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In order to visualize the mechanization of equations (B-l) - (B-3) 
for a particular application, consider the following example. 

Assume that we wish to establish optimal launch conditions from 
Earth such that minimum perihelion is achieved subject to constraints 
on the launch energy at Eartn and the transfer angle to perihelion. 
While the solution for this problem can be easily established without 
utilizing an optimization procedure, the techniques involved are repre- 
sentative of those for more complicated problems. 

The equations required _?or the assumed problem are given in 
appendix A by (A-l), (A-2), (A-3)> and (A-17). For this problem the 
variables (x's) are the velocity (V^ ) and orientation (© i ) of the 
launch trajectory. The function to be optimized (f) is the perihelion 


distance (r p ) associated with the launch trajectory. Constraints 
■^d 

include the transfer angle (q^) as well as an upper limit on launch 
energy (V^). If we assume that the desired value for transfer angle 
is 120° (equality constraint), and that V he can not exceed 10 km/sec 
(inequality constraint) then f is given by 


V 


I [°1 (T >i • 


120 ° ) + °2 ( \eo y 


where 


V = V - 10 for V > 10 
hec he he 


for V. < 10 
he — 
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and required gradient of f and partials of the constraints are given by 


Vf = 




/S s ii \ 


^he 

\6©. , dv 

\ i he, 


' 


(B-5) 


while the partial derivatives of equations (B-5) can be determined 
analytically for the simple equations involved in the assumed example, 
for large systems of equations it is usually easier to obtain numerical 
partials by perturbing the variables about some nominal value. 

Additional quantities require to initiate the iterative, optimiza- 
tion process are nominal values for the variables V, and 0., as well 

he i 

as the weighting factors or. the constraints (G's). Convergence of the 
optimization procedure is dependent on these nominal values and weighting 
factors, since unrealistic values may result in computational difficul- 
ties as well as convergence to a false optimum. For example, if the 
G's are chosen such that the contribution to f by the constraints is 
negligible as compared to unconstrained function f, then the procedure 
will tend to ignore the constraints. Likewise, if the constraints are 
weighted too heavily, then the unconstrained function f will have 
little effect on the dexermination of f . Poor choices of nominal values 
vor the variables can lead to slow convergence to the optimum resulting 
in increased computational time. 
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TABLE I.- PHYSICAL CONSTANTS 
Sun Venus Jupiter 

Earth 

u, km^/sec 2 

1.321*95 x 10 11 3.21*23 X 10 5 

1.265 X 10 8 

3.9858 X 10 5 

U, km/sec 

35.003 

15.051 

29.77 

p f km 

6.159^ X 10 5 

4.8079 x 10 T 

9.2391 X 10 5 

Distance to 
sun, km 

1.0814 X 10 8 

7.7782 x 10 8 

1.495 X 10 8 


i 
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TABLE II . - CONSTANTS USED IN EQUATIONS OF APPENDIX A 


FOR DETERMINING TRANSFER ANGLES AND TIMES 


Pre -Encounter 


Post -Encounter 


Transfer Condition Transfer Condition 


K 1 *2 K 3 K l* S 


Inside Earth 
orbit (9^ = 0) 


1 -1 


Outside Earth 
orbit (0 ± = 180°) 


-1 


-1 


Pre -Perihelion 


1 0 


Pre -Aphelion 


-1 


1 


Unperturbed 
Vehicle Orbit 


V 



Figure 1.- Illustration of encounter involving increase 

in vehicle velocity. 





Figure 2.- Illustration of various phases of typical 

flyby mission. 
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POSITIVE HELIOCENTRIC 
FLIGHT PATH ANGLE (7 g ) 
AFTER ENCOUNTER 



+ 



NEGATIVE HELIOCENTRIC 
FLIGHT PATH ANGLE (-7. ) 
AFTER ENCOUNTER 


Figure 6.- Geometry used to establish transfer times and angles for 

post -encounter phase of mission. 



Figure 7.- Geometry used to establish Earth and planet position with 
respect to post -encounter aphelion and perihelion. 
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Figure 16.- Transfer time from insertion to perihelion for various 
launch velocities and flyby distances at Jupiter. 
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Figure 17*- Earth position with respect to post -encounter perihelion for various 

launch velocities and flyby distances at Jupiter. 
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(a) Corridors for flyby distance greater than 6 planetary 

radii (fig. 1^). 

Figure 21.- Variation in approach and perigee corridors with percent 
of minimum perihelion achieved for Jupiter flyby- solar probe 
missions. 
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Figure 25.- Variation in flyby parameters with distance of closest 
approach at Jupiter for the V^ e = 10.8 km/aec condition 
yielding maximum encounter corridors (fig* 22 ). 




